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Abstract 

We present charged-particle multiplicities as a function of pseudorapidity and collision central- 
ity for the 197 Au+ 197 Au reaction at y / s/v r /v=200 GeV. For the 5% most central events we obtain 
dN c h I dr]\ v= Q = 625±55 and N c h\-4.7<ri<.4.7 = 4630 ±370, i.e. 14% and 21% increases, respectively, 
relative to y / s]v r /v=130 GeV collisions. Charged-particle production per pair of participant nucle- 
ons is found to increase from peripheral to central collisions around mid-rapidity. These results 
constrain current models of particle production at the highest RHIC energy. 

PACS numbers: 25.75.Dw 
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A central question in the study of collisions between heavy nuclei at the maximum energy 
of the RHIC facility, v /JJv r /v=200 GeV, is the role of hard scatterings between partons and 
the interactions of these partons in a high-density environment. A reduction in the number 
of hadrons at large transverse momentum has already been observed for a/stvat=130 GeV 
collisions that may hint at suppression of hadronic jets at high matter densities Jl], [|. More 
generally, it has been conjectured that new phenomena related to non-perturbative QCD 
may occur at the highest RHIC energy. For example, a saturation of the number of parton 
collisions in central nucleus-nucleus collisions could lead to a limit on the production of 
charged particles || £|, ||. 

The present Letter addresses these issues with the first comprehensive investigation of 
multiplicity distributions of emitted charged particles in relativistic collisions between 197 Au 
nuclei with v /J/v r /v=200 GeV. In particular, we have measured pseudorapidity distribu- 
tions of charged particles dN c h/dri in the range —4.7 < 77 < 4.7 as a function of colli- 
sion centrality. The pseudorapidity variable r] is related to the particle emission angle 9 
with 1] = — ln[tan(#/2)]. We find that the production of charged particles at mid-rapidity 
(77 pa 0) increases by (14±4)% for the most central collisions relative to y / J/\T/v=130 GeV 
collisions |J 0, @, [|, in agreement with results of the PHOBOS experiment |]IU|| . In highly 
energetic nuclear collisions, charged particles can be produced by hadronic ("soft") as well 
as partonic ("hard") collision processes. By extending the dN c h/dr) systematics to cover a 
range of reaction centralities and pseudorapidities, it becomes possible to more fully explore 
the different reaction mechanisms. 

The data were obtained using several subsystems of the BRAHMS experiment at 



RHIC [p j: the Multiplicity Array (MA), the Beam-Beam Counter (BBC) arrays, and the 
Zero-Degree Calorimeters (ZDCs). An analysis of charged-particle multiplicities for Au+Au 
reactions at v /J/v r /v=130 GeV using a nearly identical procedure to the one presented here is 
described in ref. |§. 

The MA determines dN c h/drj around mid-rapidity with a modestly segmented Si-strip- 
detector array (SiMA) surrounded by an outer plastic-scintillator tile array (TMA) in a 
double, hexagonal-sided barrel arrangement. Each of the 25 Si detectors (4 cm x 6 cm x 
300 /im) is located 5.3 cm from the beam axis and is subdivided along the beam direction 
into seven active strips. The TMA has 35 plastic-scintillator tiles (12 cm x 12 cm x 0.5 cm) 
located 13.9 cm from the beam axis. The effective coverage of the MA is —3.0 < i] < 3.0. 
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The SiMA is used alone for determining dN c h/dr] values near mid-rapidity because of its 
higher segmentation. However, both the SiMA and TMA are used for establishing reaction 
centrality, as discussed below. Particle multiplicities are deduced from the observed energy 
loss in the SiMA and TMA elements by using GEANT simulations [12] to relate energy loss 
to the number of particles hitting a given detector element |J. SiMA and TMA elements 
are calibrated using low-multiplicity events where well-defined peaks are observed in the 
individual energy spectra corresponding to single-particle hits [[J. 

The BBC arrays consist of two sets of Cherenkov UV-transmitting plastic radiators cou- 
pled to photomultiplier tubes. The Cherenkov radiators are positioned around the beam 
pipe with one set on either side of the nominal interaction point at a distance of 2.20 m. 
The time resolution of the BBC elements permits the determination of the interaction point 
with an accuracy of ~ 0.9 cm. Charged-particle multiplicities with 2.1 < \rj\ < 4.7 are de- 
duced from the number of particles hitting each detector, as found by dividing the measured 
detector signal by the average response of the detector to a single particle. 

The ZDCs are located ±18 m from the nominal interaction point and measure neutrons 



that are emitted at small angles with respect to the beam direction [13]. Clean selection 
of minimum-biased events required a coincidence between the two ZDC detectors and a 
minimum of 4 "hits" in the TMA. It is estimated that this selection includes 95% of the 
Au+Au total inelastic cross section. 

Reaction centrality is determined by selecting different regions in the total multiplicity 
distribution of either the MA or BBC arrays. The distributions are adjusted for "missed" 
events as described in ref. ||. In determining dN c h/drj, the centrality dependence of the MA 
and BBC distributions are based on the total multiplicity measurements of the corresponding 
array, thus allowing a range of vertex locations to be used in the BBC analysis beyond the 
acceptance of the MA (see ref. ||). For 3.0< |?7| <4.2, where it was possible to analyze the 
BBC data using both centrality selections, the two analyses give results to within 2% of each 
other. In general, statistical errors on the measurements are less than 1%, with systematic 
errors of 8% and 10% for the SiMA and BBC arrays, respectively. The systematic errors 
are dominated by overall scaling uncertainties resulting from the calibration procedures and 
should primarily lead to a common scale offset for data obtained at the two RHIC energies. 
However, there may be as much as a 3% relative scale error between the two energies. A 
point-to-point error is also present, as indicated by the small asymmetry seen in Fig. [I] for 
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FIG. 1: Distributions of dN c h/drj for centrality ranges of, top to bottom, 0-5%, 5-10%, 10-20%, 
20-30%, 30-40%, and 40-50%. The SiMA and BBC results are indicated by circles and triangles, 
respectively. Statistical errors are shown for all points where they are larger than the symbol size. 

the more central collisions. 

Figure p] shows the measured dN c h/dr] distributions for charged particles for several cen- 
trality ranges. The dN c h/dr] values for these selected centralities at r/=0, 3.0, 4.5 are listed 
in Table |, together with the average number of participating nucleons {N part ) estimated 
from the HIJING (heavy-ion jet interaction generator) model |TJJ using default parame- 
ters. For the most central collisions (0-5%), dN c fi/dr)\ v=0 = 625 ± l(stat) ± 55(syst). This 
gives a scaled multiplicity value of (dN c h/dr]) / (N part /2) = 3.5 ± 0.3 charged particles per 
participating nucleon pair and indicates a (13 ± 4)% increase relative to Au+Au reactions 
at v /J/v r /v=130 GeV |5|, |I5|]. For the most peripheral collisions analyzed here (40-50%), we 
find dN ch / drj\ v= Q = 110 ± 10, resulting in a scaled value of 3.0±0.4. By integrating the 0- 
5% multiplicity distribution we deduce that 4630 ± 370 charged particles are emitted in the 
considered pseudorapidity range. This value is (21 ±4)% higher than at ^5^=130 GeV §. 

While the scaled multiplicities increase with centrality at mid-rapidity, Fig. ^| shows 
they are independent of both collision centrality and beam energy over a pseudorapidity 
range from 0.5 to 1.5 units below the beam rapidity. This is found for energies ranging 
from the CERN-SPS energy ( v /J/v r /v=17 GeV) [16| to the present RHIC beam energy and 
is consistent with a limiting-fragmentation picture in which the excitations of the fragment 
baryons saturate at a moderate collision energy, independent of system size |J . The increased 
projectile kinetic energy is utilized for particle production below beam rapidity, as evidenced 
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TABLE I: Table of dN c h/dr] values. Total uncertainties, dominated by the systematics, are 
indicated. The average number of participants (N part ) and collisions (N co u) is given for each 
centrality class, with combined model and experimental uncertainties. The model uncertainties 
are obtained by varying the assumptions for the Glauber picture, including the NN cross section 
and the nuclear radius and surface diffuseness parameters. N c h is the integral charged-particle 
multiplicity within —4.7 < r] < 4.7. 

Cent- rj = rj = 3.0 ^ = 4.5 N ch {N coll ) {N part ) 

rality 

0-5 625±55 470±44 181 ±22 4630±370 1000±125 357±8 
5-10 501±44 397±37 156±18 3810±300 785±115 306±11 
10-20 377±33 309±28 125±14 2920±230 552±100 239±10 
20-30 257±23 216±17 90 ±10 2020±160 335±58 168±9 
30-40 174±16 149±14 64 ±7 1380±110 192±43 114±9 
40-50 110±10 95±9 43 ±5 890±70 103±31 73±8 



by the observed increase in the scaled multiplicity for central events at mid-rapidity. 

Figure |3| presents the dN ch /drj distributions obtained by averaging the values for negative 
and positive pseudorapidities to further decrease the experimental uncertainties. The solid 
lines are calculations using the model of Kharzeev and Levin ||. This model, which is 
based on a classical QCD calculation using parameters fixed to the v /saF/v=130 GeV data, 
is able to reproduce the magnitude and shape of the observed multiplicity distributions 
quite well. The dashed lines are the results of a multiphase transport model (AMPT) [17|, 
|IH| . This is a cascade model based on HIJING [fPH , but includes final-state rescattering of 
produced particles. The AMPT model is also able to account for the general trend of the 
measured distributions, particularly for the most central collisions. We also plot the similar 
distributions |19| from pp collisions at -y/i=200 GeV, scaled up by the corresponding number 
of Au±Au participant pairs, for the 0-5% and 40-50% centralities. For central collisions the 
Au+Au data show a strong enhancement over the entire pseudorapidity range relative to the 
pp results, with an excess of (41 ±9)% observed at mid-rapidity. This suggests significant 
influence of the extended, high-density medium in the case of the heavy-ion collision. The 
observed enhancement decreases to about 10% for 40-50% centrality collisions. We also 
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FIG. 2: Charged-particle multiplicities normalized to the number of participant nucleon pairs 
(see Table I) for the present 0-5% central (open circles) and 40-50% central (open squares) Au+Au 
results at ^5^=200 GeV, the BRAHMS 0-5% Au+Au results f§ at y/sjf}f=130 GeV (closed 
circles) and the 9.4% central Pb+Pb data at y / J/v r /v=17 GeV(closed triangles) of ref [16]. Data are 
plotted as a function of the pseudorapidity shifted by the relevant beam rapidity. Representative 
total uncertainties are shown for a few Au+Au points. 

note that the width in pseudorapidity of the measured distributions increases slightly as 
the centrality decreases, with a RMS = 2.33 ± 0.02 and 2.40 ± 0.02 for the 0-5% and 40-50% 
centralities, respectively. This again suggests increased particle production at mid-rapidity 
for more central collisions. These values can be compared to urms — 2.38 ± 0.05 for the pp 
data. 

The ratios of dN ch /drj values measured at v /J/v r /v=130 GeV and v /J/v r /v=200 GeV for 
different centralities are shown in Fig. An increase in charged-particle multiplicity as a 
function of energy for a central-plateau region (|?7| < 2.5) is observed, with a comparable 
increase of 10% to 20% observed for all centralities. The upturn in the ratios seen at forward 
rapidities is due to the widening of the multiplicity distribution at the higher energy, con- 
sistent with the increase in beam rapidity (Ay = 0.43). The curves show the corresponding 
ratios resulting from the two model calculations. 

Finally in Fig. |5] we plot (dN c h/drj)/{N part /2) as a function of the average number of 
participants (N part ) for three narrow pseudorapidity regions (Arj ~ 0.2) around rj =0, 3.0 and 
4.5. As already suggested, particle production per participant pair is remarkably constant 
and near unity at the forward rapidities characteristic of the fragmentation region, while 
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FIG. 3: (a-d) Measured dN c h/dr] distributions for centrality ranges of 0-5%, 5-10%, 20-30% and 
40-50%. Theoretical predictions by Kharzeev and Levin (solid line) and by the AMPT model 
(dashed line) are also shown. Result from pp collisions at -^=200 GeV p9| , scaled by the Au+Au 
values of (N part )/2, are shown with stars (a,d). 
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FIG. 4: Ratio of dN cn /dri values at y / sjv r /v=200 GeV and 130 GeV compared to the model 
calculations (see Fig. 3 caption). Total uncertainties are shown, assuming a 3% relative scaling 
uncertainty between the two energies. 



showing a significant increase for the more central collisions with r\ rs 0. The mid-rapidity 
behavior has been attributed to the onset of hard-scattering processes which are dependent 
on the number of binary nucleon collisions N co u rather than N part ||. Using N co u and 
Npart values from HIJING H] we fit the data with a function of the form dN cn /dr] = 
a ■ N part + (3 ■ N coU . For rj = 0(4.5) we obtain: a = 1.26 ± 0.09 ± 0.20 (0.66 ± 0.03 ± 0.10) 
and (3 = 0.15 =F 0.04 =p 0.05(-0.06 =F 0.01 =F 0.03), where the first uncertainty assumes a 3% 
point-to-point error for the dN cn /drj values and the second uncertainty results from the N co u 
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FIG. 5: dNch/dr] per participant nucleon pair as a function of the average number of partici- 
pants (see table) for rj= (circles), 3.0 (squares) and 4.5 (triangles). The curves show the model 



predictions (see Fig. 3 caption). The star denotes the pp result at rj = [19] 



and N part uncertainties. For comparison we find a = 1.24 ± 0.08 ± 0.20 (0.55 ± 0.02 ± 0.09) 
and = 0.12 =p 0.04 =p 0.06(-0.09 =f 0.01 T 0.03) at v /s^7=130 GeV. For central events at 
7] pd we find that the hard-scattering component to the charged-particle production remains 
almost constant, with values of (20±7)% and (25±7)% at v /^=130 GeV and 200 GeV, 
respectively. For this comparison, only the experimental component of the uncertainties are 
given since the theory uncertainties will be highly correlated at the two energies. 

In conclusion, we find that the charged-particle production increases by a constant amount 
from v /J/v r /v=130 GeV to y / J/v r /v=200 GeV in a wide region around mid-rapidity. The data 
are well reproduced by calculations based on high-density QCD and by the AMPT/HIJING 
microscopic parton model. A phenomenological analysis in terms of a superposition of 
soft- and hard-scattering particle production indicates that the hard-scattering compo- 
nent seen at mid-rapidity for central collisions is not significantly enhanced as compared 
to v /sjv r /v=130 GeV results. We find good consistency with the gluon saturation model of 
Kharzeev and Levin, but stress that within errors of models and data alike, the data can 
be equally well reproduced by other models that do not require parton-collision saturation. 
This work establishes a baseline for particle production at the maximum energy currently 
available for nucleus-nucleus collisions. 
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